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The effect of mechanical stirring on horizontal convection 
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ABSTRACT 

The theoretical analysis of the energetics of mechanically-stirred horizontal con- 
vection for a Boussinesq fluid yields the formula: 



G(APE) = lmixing G{KE) + (1 + T 



mixing 



)W r , 



laminar 



where G(APE) and G(KE) are the work rate done by the buoyancy and mechani- 
cal forcing respectively, ^mixing is the mixing efficiency, and W r ,iaminar is the back- 
ground rate of increase in gravitational potential energy due to molecular diffusion. 
The formula shows that mechanical stirring can easily induce a very strong buoyancy- 
driven overturning cell (meaning a large G(APE)) even for a relatively low mixing ef- 
ficiency, whereas this is only possible in absence of mechanical stirring if 



mixing 



> 1. 

Moreover, the buoyancy-driven overturning becomes mechanically controlled when 
JmixingG(KE) > (1 + ^1 m ixingW r , laminar ■ This result e xplains why the buoyan cy- 
driven overturning cell in the laboratory experiments by Whiteh ead et al.l (|2008f ) is 
amplified by the lateral motions of a stirring rod. The formula implies that the ther- 
modynamic efficiency of the ocean heat engine, far from being negligibly small as is 
commonly claimed, might in fact be as large as can be thanks to the stirring done by 
the wind and tides. These ideas are further illustrated by means of idealised numerical 
experiments. A non-Boussinesq extension of the above formula is also given. 



1 Introduction 



A fundamental objective of ocean circulation theory is 
to characterise and quantify the relative importance of the 
mechanical and thermodynamical forcing in driving and stir- 
ring the oceans. Although the oceans are highly nonlinear 
a system, ocean circulation theory nevertheless historically 
evolved under the assumption that separate theories for the 
wind-driven and buoyancy-driven circulation could be de- 
veloped, as it the two kind of circulation could be regarded 
as somehow decoupled. From a thermodynamic viewpoint, 
this is somehow justified from the fact that mechanical forc- 
ing, unlike the buoyancy forcing, does not modify the fluid 
parcels' buoyancy, suggesting that distinct effects should 
be observable. In the thermodynamic engineering literature, 
mechanical forcing is commonly regarded as "shaft" work, 
which alters the energy of the fluid without altering its total 
entropy. In oceanographic textbooks, wind-driven circula- 
tion theories usually occupy a more prominent place that 
buoyancy-driven circulation theories. Indeed, it turns out to 
be relatively easy to link the near surface transport (Ekman 
theory) or the vertically-integrated transport (Sverdrup the- 
ory) to the wind stress or its spatial derivatives in a quan- 
titative relatively accurate way. Theories for the buoyancy- 
driven circulation, on the other hand, are often more quali- 
tative, and at best offer only tentative scaling arguments to 
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relate the strength of the large-scale ove rturning cell taking 
place in the meridional/vertical plane, e.g. I Colin de Verdierd 
(|l993t ). From a physical viewpoint, the two main salient in- 
gredients of the buoyancy-driven circulation have been tra- 
ditionally associated with high-latitude cooling, envisioned 
as the destabilising mechanism setting up the meridional 
overturning circulation into motion, and turbulent diapyc- 
nal mixing, envisioned as the process required to carry heat 
from the surface at the depths cooled by deep water forma- 
tion, thus precluding the oceans to fill up with dense and 
salty waters. 

Physically, however, the buoyancy-driven circulation 
can only be regarded as independent of the mechanical forc- 
ing if one can establish that the stirring required to maintain 
turbulent diapycnal mixing is predominantly sustained by 
the work rate done by the surface buoyancy fluxes G(APE) 
(i.e., the rate of available potential energy production, see 
lLorenj (|l955t )h A crucial question, therefore, is how large is 
G(APE)1 In thei r observational analysis of ocean energetics, 
lOort et all (Il994n found G(APE) = 1.2 ± 0.7TW, and con- 
cluded that the buoyancy forcing was as important as the 
wind in forcing and stir ring the oce a ns. T his result, how- 
ever, was challenged by iMunk et all Jl99ct) (MW98 ther e- 
after), who contended that the so-called Sandstr5m| (|l908l )'s 
theorem forbids G(APE) to be significant, and hence that 
turbulent diapycnal mixing must be primarily be sustained 
by the work rate done by the mechanical forcing due to the 
wind and tides, leading to the idea that the buoyancy-driven 
circulation is actually mechanically-driven. MW98 further- 
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more sought to quantify the magnitude of the mechanical 
sources of stirring required to sustain diapycnal mixing in 
the oceans. Their main result is the following formula: 



G(KE) = 



Wr, forcing 

5 

^mixing 



(i) 



where G(KE) is the work rate done by the mechanical 
sources of stirring, W r , forcing is the rate at which the back- 
ground gravitational potential energy GPE r decreases as 
the result of high-latitude co oling, and '"/mi xing is the so- 
called mixing efficiency, e.g., lOsbornl l|l980l ). MW98 esti- 
mated W r j rcing ~ 0.4 TW from assuming the deep water 
formation rate to be known, and used the canonical value 
^mixing = 0.2 for the mixing efficiency to conclude that 
G(KE) — 0(2TW) was required to sustain the observed 
rate of turbulent diapycnal mixing and its associated pole- 
ward heat transport of about 2 PW. Since the wind forcing 
provides no more than about 1 TW, this result suggests an 
apparent shortfall of about 1 TW of mechanical stirring that 
MW98 argued must only come from the work rate done by 
the tides, spawning much research over the past decade on 
the issue of tidal mixing. 

MW98's study prompted much debate over the past 
decade about the relative importance of the surface buoy- 
ancy forcing, and about whether Sandstrom's the orem re- 
ally im plied for G(APE) to be negligible. Recently. ITailleuxl 
(2009) revisited a number of misconceptions about the 
nature of energy conversions in turbulent stratified fluids 
by extending the av ailable potential energy framework of 
IWinters~et al. ( 1995) to the fully compressible Navier-Stokes 
equations (CNSE thereafter). One of the main result was to 
demonstrate from first principles that: 



G(APE) « W r , 



forcing j 



(2) 



which physically states that the production rate of APE is 
approximately equal to the rate at which GPE r decreases as 
the result of high-latitude cooling. This result is extremely 
important, because it rigourously proves that it is inconsis- 
tent to assume W r j rcmg to be large while simultaneously 
assuming G(APE) to be negligible, as done in MW98. As 
a result, rather than establishing that the work rate done 
by surface buoyancy forcing is small, MW98's study actu- 
ally establishes precisely the contrary, since their estimate 
Wrjorcing ~ 0.4 TW impl ies G(APE) » 0.4 TW, which 
turns out to be very close to | Oort et aL (1994)'s lower bound 
for G(APE). Moreover, ITailleuxl i|2009l ) was able to show 
that a generalisation of MW98's Eq. {1} for a non-Boussinesq 
fluid can be rigourously derived from first principles by com- 
bining the mechanical energy budget for APE and GPE r 
separately, the final result taking the form: 



G{KE) = 



1 + (1 - 07" 



C7n 



-G(APE) 



G(APE), 



(3) 



where £ is a non-Boussinesq nonlinearity parameter which is 



such that — oo < £ ^ 1 for water and seawater, j„ 



is the 



mixing efficiency, and Rf = "/mixing I + "/mixing) is the so- 
called flux Richardson number, e.g.. lOsborn1 (| 19801 ). MW98's 
Eq. fTJ is recovered in the limiting case f = 1 describing a 
Boussinesq fluid with a linear equation of state, by using the 
above result that G(APE) « W r , forcing- 



The fact that the generalisation of MW98's result can 
be written as a formula linking G(KE) and G(APE) rather 
than G(KE) and W r , forcing raises the question of whether it 
is really appropriate to interpret Eq. © (and hence Eq. {T}) 
as a constraint on the amount of G(KE) required to sustain 
diapycnal mixing in the oceans, as suggested by MW98. In- 
deed, it seems that such an interpretation implicitly assumes 
that both G(APE) and "/mixing can be regarded as fixed in 
some sense, but whether this is the correct way to interpret 
Eq. ([3]) is unclear, given that the latter can be rewritten in 
the two following ways: 



G(APE) = 



£7" 



1 + (1 - 0~fn 



-G(KE), 



G(APE) 



(4) 



(5) 



G(KE) + G(APE) ' 

i.e., under the form of a constraint either on G(APE) or 
on £,Rf. Note that from observations, plausible numbers for 
the wind and buoyancy forcing are G(KE) ~ 1 TW and 
G(APE) « 0.5 TW, which, if inserted into Eq. © yields: 



05 
1.5 



= 0.3333. 



Another way to pose the problem initially formulated by 
MW98 would be to ask the question of whether the above 
value for £Rf can be ruled out from what we know about 
mixing efficiency (and the role played by the nonlinearity of 
the equation of state, but we know less about the latter). 
Although the value "/mixing ~ Rf ~ 0.2 is often assumed, 
it is important to realize that such value mostly pertains to 
mechanically-driven turbulent mixing, as occurs in relation 
with shear flow instability for instance. Indeed, mixing effi- 
ciency valu es as high as Rf = 0.5 ("/mixing = 1) have been 
reported bv lDalziel et al.l (|2008l ) for buoyancy-driven mixing 
associated with Ray leigh- Taylor instability. Since G(APE) 



is comparable to G(KE) , values of "/„ 



intermediate be- 



tween those for mechanically- and buoyancy-driven mixing 
should be expected. 

In trying to figure out the exact physical meaning of Eq. 
((3]), it is important to discuss the nature and structure of 
G(KE) and G(APE) in more details. Thus, the work rate 
done by the wind stress is given by: 

G(KE) = t- u s dS (6) 
Js 

where r is the wind stress at the ocean surface, and u s is the 
ocean surface velocity. The important point about G(KE) is 
that it is a correlation between the external forcing (the wind 
stress) and a parameter depending upon the particular state 
of the system (the surface velocity). In other words, it is es- 
sential to realize that G(KE) cannot be determined from the 
knowledge of the external forcing alone. This is important, 
because Eq. (JS} expresses the possibility for the wind-driven 
circulation to be controlled by the buoyancy forcing, to the 
extent that the latter is able to influence the surface veloc- 
ity. In order to gain insight into the potential importance 
of the buoyancy control of the wind-driven circulation, it 
would be interesting to compute G(KE) for a purely wind 
forced homogeneous ocean model. Doing so, however, is be- 
yond the scope of this paper. For the present purposes, we 
shall assume that G(KE) is primarily determined by the 



MECHANICALLY-STIRRED BUOYANCY-DRIVEN CIRCULATIONS 3 



wind forcing, and that buoyancy forcing only alters it as a 
second order effect, which seems plausible on the basis that 
the near surface circulation in the oceans is usually thought 
to be the signature of the wind forcing rather than the buoy- 
ancy forcing. 

With regard to the work rate done by surface buoyancy 
fluxes, it was shown by ITailleuxl (|2009h to be given by the 
following expression: 



G{APE) = 



T-T r 
T 



'surf 



dS, 



(7) 



in the case of a compressible thermally-stratified fluid forced 
by surface heat fluxes, where T is the surface temperature of 
the fluid parcels, and T r the temperature the surface parcels 
would have if displaced adiabatically to their level in lLorena 
l|l955l Vs reference state, while Q 3ur f is the diabatic rate of 
heating cooling/heating due to the surface heat fluxes. Use- 
ful approximations for G(APE) can be obtained by expand- 
ing T as a Taylor series expansion around the surface pres- 
sure P a , i.e., T^T r + F r {P a - P r ) + 0((P a - Pr) 2 ), where 
T r = a r Tr/(p r Cpr) is the adiabatic lapse rate, leading to: 



G(APE) 



a r {P r - Pa 

prCpr 



{ surf 



dS 



a r gz r 

Cpr 



Q sur fdS , 



(8) 



by using the approximation P r — P a ~ —pogz r , the lat- 
ter expression being actual ly the Boussinesq a pproxima- 
tion of G(APE) derived by IWinters et"aL1 (|l995h provided 
that one regards the thermal expansion a r and specific 
heat capacity C pr as constant (the suffix r means that the 
variables have to be estimated in lLoren i (|!955lVs refer- 
ence state). As for G(KE), G(APE) also takes the form 
of a correlation between the external forcing (the heat- 
ing/cooling rate at the surface) and a parameter depending 
on the particular state of the system, namely (T — T r )/T ~ 
&r(Pa Pr) / (PrGpr ) » a r gz r /C P r- Eqs. (J7J) and flSJ) make 
it possible, therefore, for the buoyancy-driven circulation to 
be mechanically-controlled. Physically, this is expected, be- 
cause mechanical forcing is widely agreed to increase diapy- 
cnal mixing, which should produce a deeper thermocline, 
thereby allowing dense plumes to penetrate deeper, thus 
increasing G(APE) and hence the buoyancy-driven circu- 
lation. Interestingly, th is appears consistent wi th the labo- 
ratory experiments by [Whitehea d et af] (|2008T ) which pro- 
vide evidence that the lateral motion of a stirring rod can 
greatly enhance the strength of horizontal convection (see 
iHughes et ail (|2008l l for a review on this topic). This paper 
examines the possibility of interpreting this result as the 
consequence that G(APE), and hence the buoyancy-driven 
circulation, is increased by the work rate done by the stir- 
ring rod. To that end, Section 2 seeks to extend Eq. (J3J) to 
also describe the purely buoyancy-driven case for which it is 
not valid. The physical meaning of the extended formula is 
then explored by means of idealised numerical experiments 
in Section 3. Finally, section 4 presents a discussion of the 
results. 



2 Theoretical results about the effects of 

mechanical stirring on the work rate done by 
surface buoyancy fluxes 

2.1 Summary of Tailleux (2009) 's theory 

Our starting point is the theoretical description of the 
energetics of mechanically and thermodyn amically forced 
turbulent stratified fluids recently derived by ITailleuxl (|2009T l 
pertaining to a thermally-stratified fluid governed by the 
fully compressible Navier-Stokes equations (NCSE there- 
after), which builds upon the avail able potential energy 
framework previously introduced by IWinters et al.l (|l995T l 
for a Boussine sq fluid with a linear equation of state. In 
iTailleuxl (|2009T ) 's description, the energetics is described at 
leading order by means of the following five evolution equa- 
tions: 



d(KE) 

~dt " 
d(APE) 
It 

d(GPE r 
dt 

d(IEp) _ 
dt 

d(^I Eexer 

~dT~ 



-- -C{KE, APE) + G(KE) - D(KE) , (9) 
= C(KE,APE) + G{APE) - D{APE), (10) 



r .laminar 



W, 



r. forcing i 



(11) 



(l-r )Q net +D(APE)+D(KE)-G(APE), (12) 



= ToQr, 



w, 



r .turbulent 



Wlar, 



(13) 



w here KE is t he volume-integrated kinetic energy, APE 
is lLorenzl l|l955h 's volume-integrated available potential en- 
ergy, GPE r is the vo l ume-in tegrated gravitational poten- 
tial energy of iLorenzl ljl955i )'s reference-state, IEo is the 
volume-integrated of a subcomponent of internal energy 
(IE) that we call the dead part of IE, IE exergy is the 
volume-integrated of another subcomponent of IE called the 
exergy. Physically, variations in IEo are associated with vari- 
ations in the equivalent thermodynamic equilibrium temper- 
ature To(t) of the system, whereas variations in IE exsrgy re- 
flect variations in the reference temperature profile T r (z,i). 
The other important conversion terms are C(KE , APE), 
the so-called buoyancy flux that represents the reversible 
conversion between KE and APE; G(APE), the production 
rate of APE that physically represents a conversion between 
IE a and APE; W r ,iaminar and W r ,turbuie.nt represent the 
laminar and turbulent rate of exchange between IE exergy 
and GPE r due to molecular diffusion; D(KE) is the dis- 
sipation rate of KE into internal energy IEo by molecular 
viscous processes; D(APE) is the dissipation rate of APE 
into internal energy IEo by molecular diffusive processes. 
The parameter To ^ 1 plays the role of a thermodynamic 
efficiency that is very small for a nearly incompressible fluid, 
and which controls how much of the net surface heating rate 

Qnet splits between IE exerq y and IEp. 

As previously shown bv lwinters et al-H 19951 ) . a number 
of conversion terms appear to be strongly correlated to each 
other. This is the case for D(APE) and W r ,turbuient, owing 
to both terms: 1) being controlled by molecular diffusion; 
2) being co ntrolled by the s pectra l dis tribution of th e APE 
density, see iHollidav et~aH (| 198 if ) and iRoullet et ail (|2009h 
for a discussion of the l atter concept. F or a compressible 
thermally-stratified fluid, ITailleuxl (2009) showed that this 
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correlation between D(APE) and W r ,turbuient can be writ- 
ten under the general form: 



ir, 



r .turbulent 



: iD(APE) 



(14) 



where £ is a parameter that measures the importance of 
the nonlinearity of the equation of state, such that f = 1 
for a Boussinesq fluid with a linear equation of state, but 
— oo < £ ^ 1 for water or seawater. The second type of con- 
version terms that are correlated to each other are G(APE) 
and Wrjordng owing to b oth terms being controlled by the 
surface buoyancy forcing. ITailleuxl l|2009T l shows that to a 
good approximation, one usually has: 



G(APE) » W r , forclng . (15) 

For th e sake of brevity, the reader is referred to ITailleuxl 
(2009) for the details about the explicit forms of all the 
terms entering the above energy equations, as they are unim- 
portant for the arguments developed in this paper. Fig. {TJ 
schematically illustrates the energetics of a mechanically and 
thermodynamically forced stratified fluids associated with 
the above equations. For all practical purposes, G(APE) 
and G(KE) represent the work rate done by the buoyancy 
and mechanical forcing respectively, whereas D(APE) and 
D(KE) represent the two terms dissipating the "available" 
mechanical energy ME = APE + KE. 

2.2 A general theory for G(APE) 

A fundamental question posed by this paper is what 
controls the magnitude of the buoyancy-driven circulation, 
that is, the circulation drawing its energy from G(APE), 
with and without mechanical stirring acting on the fluid. In 
order to answer this question, we first need to introduce a 
couple of parameters that are traditionally used to measure 
the efficiency of turbulent mixing, the so-called mixing effi- 
ciency ^/mixing and the flux Richardson number Rf, which 
are defined here as follows: 



"l mixing 



D(APE) 
D{KE) ' 



Rf 



"^mixing 
1+7 ■m i x t n 



D(APE) 



D(APE) + D(KE) ' 



(16) 



(17) 



The physical rationale for such definitions, as well as their 
connec tion to other ex isting definitions, was discussed in de- 
tails in ITailleuxl (2009) to which the reader is referred to for 
details. From a practical viewpoint, the differences in ex- 
isting definitions is unimportant, as least in the context of 
a Boussinesq fluid, as then all definitions are then numeri- 
cally equivalent, even if based on different physical assump- 
tions. In particular, the present definitions allow to recover 
MW98's Eq. ([1]), as disc ussed be l ow. T o distinguish them 
from existing definitions, ITailleuxl { 2009) refer to the above 
"/mixing and Rf &s the " dissipative" mixing efficiency and 
flux Richardson number respectively. 

In a second step, the mechanical energy balance is con- 
structed by summing the steady-state version of the KE 
and APE equations, which yields: 



G(APE) + G(KE) = D(APE) + D(KE). 



(18) 



Now, by combining this equation with the definition of the 
flux Richardson number (i.e., Eq. (|17[) . one may write: 



D(APE) = Rf [G(APE) + G(KE)]. (19) 
Next, we turn to the steady-state GPE r balance, viz., 



^Vr, turbulent ~t~ V^V, laminar — Vt^r, f orcing ■ 



(20) 



As mentioned above, one has to a very good approximation 
the following equality W r , forcing ~ G(APE). For mathe- 
matical rigour, we introduce a parameter £2 ~ 1 making 
exact the equality W r , forcing = £,2G(APE). Using the defi- 
nition of £ above, the GPE r budget provides the following 
relationship: 



£D(APE) + WrJaminar = &G(APE), 



which after some algebra eventually yields: 
G(APE) = 



^^G(KE) + WrM 



m 1 n a r 



6 - ZRf 



(21) 



(22) 



Eq. (I22|l is one of the most important result of this paper, 
as it represents one way to demonstrate the interconnection 
of the work rate done by the wind and buoyancy forcing in 
general. A useful limit is the case of the widely used Boussi- 
nesq model for which £ = £2 = 1, in which case Eq. 
simplifies to: 



G(APE) 



Rf 



1-R 



= 7rr 



g G(KE) + (!+"/„ 



-G(KE) + 

)W r ,lan 



W r . 



1-R f 



(23) 



Eq. (1 2 2 tl shows that in absence of mechanical forcing, the 
work rate done by the surface buoyancy fluxes is given by: 



si, Arm W r .laminar 

G{APE)= __ 

or, for a Boussinesq fluid: 



(24) 



G(APE) = Wr.laminar _ ^ + J mixing ) Wr.laminar ■ (25) 
1 - Rf 

The latter two results are interesting, because they show 
that although W r ,iaminar is small by construction, it does 
not forbid in principle G(APE) to be large provided that 
&—£Rf or 1 — Rf can become small enough. The problem, 
however, is that the latter conditions require Rf to be close 
to unity for a Boussinesq fluid, or equivalently "/mixing 2> 1, 
which is much larger than the widely used canonical value 
"/mixing ~ 0.2. On the other hand, we are not aware of any 
published value for "/mixing pertaining to horizontal convec- 
tion, so that sufficient empirical evidence to be really conclu- 
sive. The idealised numerical experiments presented below 
seeks to get insight into this issue. 

Perhaps the most important feature of Eq. (|22[) is to re- 
veal that even a small amount of mechanical forcing is suffi- 
cient to radically alter the nature of horizontal convection, as 
it suggests that the latter becomes mechanically controlled 
when G(KE) is such that ^RfG(KE) S> W r ,iaminar, which 
is possible even for relatively low values of "/mixing — a very 
important point. This theoretical result suggests therefore 
that G(APE) can be dramatically enhanced by the presence 
of mechanical stirring, provided that the latter provides pos- 
itive work to the system (i.e., such that G(KE) > 0). Neg- 
ative work, on the other hand, should reduce the strength 
of the buoyancy-driven circulation. The numerical examples 
studied next seek to illustrate these different ideas. 
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KE 



IEo 
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C(KE,APE) 
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r, forcing 



w, 



r,mixing 



GPEr 



Figure 1. Schematic diagram of the energy conversions taking place in a mechanically and thcrmodynamically forced turbulent thermally- 
stratified fluid corresponding to Eqs. (19113ft previously derived by Tailleux (2009). 




3 Numerical experiments 

The physical implications of the formula given by Eqs. 
(|22[) and (|23[) are explored in the following by means of nu- 
merical experiments aiming at illustrating the main effects 
mechanical stirring can have on buoyancy-driven circula- 
tions. 

3.1 Model description 

To that end, the idealised p roblem of horiz o ntal c on- 
vection previously considered bv lPaparella et ail (I 2002h is 
investigated here, the novelty being the addition of mechan- 



ical stirring to the problem modelled here as a forcing term 
in the vorticity equation. The numerical i mplementation o f 
such a model is that previously described bv lMarchall (|2007T ) ■ 
The equations solved by the numerical model are the Boussi- 
nesq equations for a fluid with a linear equation of state: 

^ = J(*,Lj)-ga^+»V 2 u J + F(x,z,t) (26) 
— = J(*,T) + K V 2 T, (27) 

where lu = V 2 ^ is the vorticity, T is the temperature, 
J(a,b) — da/dxdb/dz — da/dzdb/dx is the Jacobian op- 
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Figure 3. Temperature field for the different experiments: (a) No mechanical forcing; (b) Clockwise mechanical forcing; (c) Weak 
anti-clockwise mechanical forcing; (d) Strong anti-clockwise mechanical forcing. 



erator, g is the acceleration due to gravity, a is the thermal 
expansion coefficient assumed to be constant, v is the kine- 
matic viscosity, k is the thermal diffusivity, and F is a term 
aimed at mod elling the effect of mechanical stirring. 

Following Marchal (2007j), the above system is made 
dimensionless as follows: t = (L 2 / ' n)t* , (x, z) = L(x,,z,), 
to = (k/L 2 )uj„, T = ATT,, and * = k*„, where the starred 
quantities are the dimensionless ones 
forms of the above equations become, 
start for clarity: 



The dimensionless 
after dropping the 



dT 
~dt 



j(*,t) + v 2 r, 



(28) 
(29) 



where R a — gaAT 3 /(vk) is the Rayleigh number, and P r = 
v/k is the Prandtl number. In the numerical experiments 
described here, we used P r = 10 and R a = 10 . For com- 
pari son, note that typical o ceanic values are R a = O(10 20 ), 

e.g., IPaparella et all (I 20021 ) 

As in IPaparella et all (| 20021 ). the fluid is forced by a 
surface temperature boundary condition varying linearly in 
x. The numerical resolution used for the experiments is 51 x 
51 in a 2-D square geometry. At equilibrium, the forcing 
term in the vorticity equation is associated with the work 
rate: 



G(KE) 



J 1 * dxdz 



(30) 



which can be in principle either positive or negative. 



3.2 Experiments 

Four idealised experiments were considered: 

(i) Purely buoyancy-driven (i.e., no mechanical forcing); 

(ii) Clockwise Mechanical forcing; 

(iii) Weak anti-clockwise mechanical forcing; 

(iv) Strong anti-clockwise mechanical forcing. 

The purely buoyancy-driven case is associated with a 
clock-wise thermally-direct circulation, and represents the 
"control" reference simulation to be compared with the 
mechanically-stirred ones. The purpose of Experiment (ii) 
is to gain insight into the case for which the mechanical stir- 
ring tends to re-enforce the existing thermally-direct clock- 
wise circulation, implying G(KE) > 0. From Eq. (I23p . the 
expectation is that the overturning circulation should in- 
crease both as the result of increased G(APE), as well as 
from the direct effect of the mechanical forcing acting as a 
source of clockwise vorticity. In such a case, therefore, KE, 
APE, *, D(APE) and G(APE) all should increase. 

In experiments (iii) and (iv), however, the anti- 
clockwise forcing works against the existing buoyancy-driven 
circulation. If too weak, as designed to be the case for (iii), 
the mechanical forcing is unable to reverse the sense of the 
existing thermally-direct circulation anywhere, resulting in 
G(KE) < 0. According to Eq. (|23|l . a decrease in all quan- 
tities KE, APE, *, D(APE), and G(APE) is expected in 
that case. If the clockwise forcing is strong enough, however, 
it can locally generate a anti-clockwise circulation associ- 
ated with G(KE) > 0, making it possible in that case for 
all above quantities to be increased. 
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G(KE) 


G(APE) 


D(APE) 


^mixing 


(a) 





5.98 


4.12 


1.78 


(b) 


7.54 


10.1 


9.52 


1.04 


(c) 


-0.77 


4.72 


2.58 


1.49 


(d) 


9.83 


8.50 


7.18 


0.60 



Table 1. Values of the work rate done by the mechanical forcing 
G{KE), rate of APE production G(APE), diffusive APE dissi- 
pation rate D(APE), and mixing efficiency "(mixing for (a) No 
mechanical forcing; (b) Clockwise mechanical forcing; (c) Weak 
anti-clockwise mechanical forcing; (d) Strong anti-clockwise me- 
chanical forcing. 





KE 


APE 


* max 


H max 


(a) 


0.13 


4.04 


7.2 


3.15 


(b) 


0.64 


6.25 


15.5 


5.09 


(c) 


0.09 


5.49 


6.08 


2.19 


(d) 


0.51 


11.1 


9.67 


2.75 



Table 2. Values of the volume-integrated kinetic energy KE, 
volume-integral available potential energy APE, maximum value 
of the thermally-direct overturning streamfunction, and maxi- 
mum value of the heat transport for: (a) no mechanical forcing; 
(b) Clockwise mechanical forcing; (c) Weak anti-clockwise me- 
chanical forcing; (d) Strong anti-clockwise mechanical forcing. 

3. 3 Results 

The results are illustrated by the figs. ((2)1 and <[3j , which 
show the streamfunction and temperature distributions for 
the 4 different experiments considered. In addition, Tables 1 
and 2 summarise how several important quantities vary as a 
function of the mechanical forcing imposed. As expected, all 
the plots for ^ show evidence of a well-marked thermally- 
direct cell in all cases. Fig. ((2| (b) shows that the strength 
of the overturning cell is greatly enhanced by the addition of 
clockwise mechanical forcing. This can be explained by the 
observed increase in G(APE) (Table 1), as well as by the di- 
rect effect of the mechanical forcing as a source of clockwise 
vorticity. Fig. [3] (b) shows that the increase in G(APE) is 
associated with the dense plumes penetrating deeper in that 
case. Fig. @ (c) shows that a weak anti-clockwise mechani- 
cal forcing decreases the buoyancy-driven circulation, which 
is expected from the theory when G(KE) < as is the case 
here (Table 1). The weakening of the circulation results from 
the decrease in G(APE) , despite an apparent increase in the 
thermocline depth (Fig. [3} which one would normally asso- 
ciate with an increase in G(APE). As shown in Table 2, this 
occurs because of a reduction in the strength of the surface 
heat flux, which results in a decrease in the maximum heat 
transport H max . Fig. ((2} (d) shows that the anti-clockwise 
mechanical forcing is strong enough to generate a thermally- 
indirect circulation, resulting in G(KE) > 0, which increases 
the overturning strength. In that case, there is no ambiguity 
that the increase in the latter is entirely explained by the 
increase in G(APE) , given that the mechanical forcing acts 
as a source of anti-clockwise vorticity, which can only tend 
to weaken the buoyancy-driven thermally-direct cell. 

In the literature, there is a tendency to regard the mix- 
ing efficiency "/mixing as some kind of universal parameter 
whose value is relatively constant and uniform throughout 
the oceans. It is interesting, therefore, to examine whether 



this is the case in the present numerical experiments. To 
that end, the viscous dissipation D(KE) and diffusive dis- 
sipation D(APE) were estimated for all experiments, the 
values being reported in Table 1. As it turns out, "/mixing 
is found to vary widely across the different experiments. Its 
maximum value is achieved in the purely buoyancy-driven 
case, for which "/mixing ~ 1.78. Interestingly, this value is 
significantly l arger than the value "/mixing ~ 1 reported by 
iDalziel et al.l ([2008) in the context of Rayleigh- Taylor in- 
stability, and which represents so far the highest value of 
mixing efficiency ever reported. The Rayleigh- Taylor insta- 
bility, however, is known to have the peculiar property that 
the APE driving the instability i s not entirely av ailable for 
turbulent mixing. For this reason, iTailleuxl (|2009h suggested 
that higher values of mixing efficiency could in principle be 
achieved for buoyancy-driven turbulent mixing events if the 
above constraint could be relaxed. As it turns out, hori- 
zontal convection does not suffer from the limitations at- 
tached to the Rayleigh- Taylor instability, making it possible 
to achieve values of "/mixing significantly large than unity. 
The addition of mechanical stirring, however, is seen in Ta- 
ble 1 to systematically reduce "/mixing , the largest reduction 
occurring for the strong anti-clockwise mechanical forcing. 
In that case, a value of "/mixing = 0.6 is reached, which is 
much closer to the value "/mixing = 0.2 commonly used in 
the literature. The present numerical experiments therefore 
suggest that the mixing efficiency of a fluid forced both me- 
chanically and thermodynamically is likely to be potentially 
highly sensitive to the particular configuration considered. 
This challenges the notion that "/mixing = 0.2 should be used 
in the oceans, given that G(APE) and G(KE) are compa- 
rable in magnitude in reality. 

Another interesting question is what is the effect of the 
mechanical stirring on the heat transport, which was also 
raised in MW98's study. The last column of Table 2 provides 
insight into this issue by providing the maximum value of the 
heat transport for each experiment. The most striking result 
is that the heat transport appears to be drastically increased 
when the mechanical forcing supports the thermally-direct 
circulation, but decreased when the mechanical forcing op- 
poses the latter. In the actual oceans, the wind forcing cre- 
ates surface Ekman cells that are alternatively helping or 
opposing the large-scale thermally direct cell in the Atlantic 
ocean responsible for the Atlantic heat transport. It seems 
therefore difficult to conclude at this stage whether the over- 
all net effect of the mechanical forcing in the oceans is to 
directly contribute to the strength of the overturning circu- 
lation by being a net source of vorticity of the right sign, or if 
its contribution is only indirect and limited to its increase of 
G(APE) and hence of the strength of the buoyancy-driven 
component of the overturning. 



4 Discussion and conclusions 

In this paper, the effect of mechanical stirring on 
buoyancy-driven flows was investigated both theoretically 
and numerically. The most important theoretical result are 
the formula Eqs. (|22[) and ([23} linking the work done by 
the mechanical stirring to the work rate done by the surface 
buoyancy fluxes via the bulk mixing efficiency of the system, 
valid for the compressible Navier-Stokes equations and for 
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the Boussinesq model respectively. T hese formula ar e a fur- 
ther generalisa tion o f that d erived bv lTailleuxl (|2009h . which 
extends iMunk et al.l (|l998h 's previous result on the energy 
requirement on the mechanical sources of stirring to sustain 
diapycnal mixing in the oceans. The actual physical mean- 
ing of the formula Eqs. (|22[l and (|23p is not entirely clear, 
however. Physically, the formula only seem to express the 
existence of a mechanical control on buoyancy-driven flows 
and vice-versa, rather than a true constraint o n G(KE), in 
contra st with the interpretation put forward bv lMunk et al.l 
(1998). This idea is further re-enforced by the structure of 
G(KE) and G(APE), which demonstrate that the work rate 
done by the mechanical and buoyancy forcing does not sim- 
ply depend on the external forcing (i.e., the wind and surface 
buoyancy fluxes) , but also on the actual state of the system. 
In other words, the work rate done by either the wind and 
buoyancy forcing depends sensitively on the work rate done 
by the other kind of forcing. In this paper, this idea is exam- 
ined from the viewpoint of the buoyancy-driven circulation, 
by demonstrating that the presence of mechanical stirring 
can have dramatic effects on the value of G(APE) owing 
to its direct effect on diapycnal mixing and hence on the 
structure of the thermocline. The other way that mechanical 
forcing can affect G(APE) is by modifying the near surface 
temperature, thereby modifying the net flux into the sys- 
tem, since the flux is also to be determined as part of the 
solution when a surface temperature boundary condition is 
imposed. 

In the purely buoyancy-driven case, Eq. (I22|l demon- 
strates that large values of G(APE) are in principle 
achievable provided that the "dissipative " mixing efficien cy 
Jmixing = D(APE) /D(KE) defined in iTailleuxj l|2009l ) is 
large enough. The particular numerical experiment consid- 
ered in this paper suggests that this is possible, since the 
value "/mixing ~ 2 was reached in the purely buoyancy-driven 
case, which is one order of magnitude larger than the value 
"/mixing = 0.2 currently assumed in most studies of turbulent 
mixing. Although the result is only tentative, since a more 
systematic investigation is required to ascertain the robust- 
ness of the results, it is important to note t hat this appears 
to be c onsistent with the numerical results of lPaparella et all 
( 2002) that show entropy production in horizontal convec- 
tion to increase with the Rayleigh number. According to our 
theoretical results, this is indeed only possible if D(APE) 
and hence G(APE) also i ncrease with th e Rayle igh num- 
ber. On the other hand, iPaparella et~ai] l| 20021 ) 's ' 'anti- 
turbulence theorem" states that the viscous dissipation of 
kinetic energy is bounded from above by a bound that is 
independent of the Rayleigh number. If D(APE) increases 
with R a with no corresponding increase in D(KE), then 
one can infer that the dissipative mixing efficiency "/mixing 
must also increase with the Rayleigh number. Testing this 
hypothesis further will be an important future research goal. 

As shown by Paparell a et al.l |j2002), the "anti- 
turbulence theorem" places a stringent constraint on the 
magnitude of the viscous dissipation rate of kinetic en- 
ergy that can be sustained by buoyancy forcing alone. As 
a consequence, mechanical forcing appears critical to pro- 
duce significant amount of viscous dissipation of kinetic en- 
ergy D(KE). However, because of the relatively low mix- 
ing efficiency of mechanically-driven mixing, it is usually 
the case that the addition of mechanical forcing to hori- 



zontal convection results in relatively modest increase in 
D(APE) compared to that in D(KE). The overall conse- 
quence is to decrease the mixing efficiency "/mixing of the 
system as compared to the purely buoyancy-driven 
verified in the numerical experiments considered here for 
which G(KE) > 0. The decrease in mixing efficiency result- 
ing from mechanically stirring horizontal convection does 
not preclude an increase in G(APE) however, as expected 
from the theoretical formula, and verified in the numeri- 
cal experiments. In our opinion, this is the fundamental 
reason why the s treng th of the overturning circulation in 
I Whitehead et af1l|2008h 's laboratory experiments appears to 
be greatly enhanced by the action of the stirring rod. De- 
pending on the particular circumstances, it appears possible 
for the mechanical forcing to contribute directly to the ob- 
served overturning increase, in addition to the indirect effect 
it has on the increase of G(APE), as in the particular case 
of the clockwise mechanical forcing experiment. In the cases 
where the mechanical forcing acts anti-clockwise, however, 
only the increase in G(APE) appears to be directly respon- 
sible for the increase in the overt urning. The latte r case is 
likely to be the one pertaining to lWhitehead et al. | (l2008l )'s 
laboratory experiments, as it seems difficult to see how the 
lateral motions of the stirring rod could contribute to the 
vorticity source required to drive a thermally direct cell. 

The present results make it hard to regard Eq. (|22[) as 
a constra int on the mechan ical sources of stirring, as pro- 
posed by IMunk et~ail (|l998h . given that G(APE) is likely 
more sensitive to a change in the mechanical forcing than 
G(KE) to a change in the buoyancy forcing, although this 
has not been directly verified. The other main difficulty in 
interpreting Eq. (|22p as a constraint on G(KE) comes from 
the fact that the value of "/mixing , far from having the value 
of 0.2 commonly affixed to it, appears to be extremely sen- 
sitive to the details of the mechanical and buoyancy forcing. 
With the present experiments, values were found to range 
between 0.6 and 1.78. The fact that buoyancy forcing is as 
important as the mechanical stirring due to the wind and 
tides has important consequences for the value of mixing effi- 
ciency one sho uld use in the ocean s. We note that if one uses 
"/mixing = 1 in IMunk et al.l l|l998h 's paper, the requirement 
on the mechanical sources of stirring reduces drastically to 
G(KE) — 0.4 TW, which is easily achieved by the wind 
alone. Future research, therefore, should aim at understand- 
ing the physical mechanisms controlling the value of mixing 
efficiency in mechanically- and thermodynamically forced 
stratified fluids. We note that to date, mixing efficiency is 
traditionally studied in the context of freely decaying tur- 
bulence, no t in the context of forc ed /dissipated systems. In- 
terestingly, |Pe Boer et al.l l|2008l ) suggest that coarse reso- 
lution numerical ocean models behave as predicted by the 
present theory, as they found no overturning circulation cell 
in the total absence of wind forcing in the Atlantic ocean. 
Another issue of interest are the consequences of mechan- 
ical control for the multiple equi libria of the thermoh aline 
circulation, rece ntly discussed by I Johnson et all (|2009l ) and 
iNof et aT] (|2007l '). Indeed, such a problem is usually studied 
under the assumption that the turbulent mixing parameters 
can be somehow regarded as fixed. The present results sug- 
gest, however, that the overall value of "/mixing might differ 
for different equilibria. 
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